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Abstract 



Realistic effective supergravity models have a variety of sources of lepton flavour vi- 
olation (LFV) which can drastically affect the predictions relative to the scenarios 
usually considered in the literature based on minimal supergravity and the super- 
symmetric see-saw mechanism. We catalogue the additional sources of LFV which 
occur in realistic supergravity models including the effect of D-terms arising from 
an Abelian U(l) family symmetry, non-aligned trilinear contributions from scalar F- 
terms, as well as non-minimal supergravity contributions and the effect of different 
Yukawa textures. In order to quantify these effects, we investigate a string inspired 
effective supergravity model arising from intersecting D-branes supplemented by an 
additional U(l) family symmetry. In such theories the magnitude of the D-terms is 
predicted, and we calculate the branching ratios for \i — > e-f and r — > for different 
benchmark points designed to isolate the different non-minimal contributions. We 
find that the D-term contributions are generally dangerously large, but in certain 
cases such contributions can lead to a dramatic suppression of LFV rates, for ex- 
ample by cancelling the effect of the see-saw induced LFV in r — > /ry models with 
lop-sided textures. In the class of string models considered here we find the surprising 
result that the D-terms can sometimes serve to restore universality in the effective 
non-minimal supergravity theory. 
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1 Introduction 



Lepton flavour violation (LFV) is a sensitive probe of new physics in supersymmetric 
(SUSY) models [HE]- In SUSY models, LFV arises due to off-diagonal elements 
in the slepton mass matrices in the 'super-CKM' (SCKM) basis, in which the quark 
and lepton mass matrices are diagonal [3]. In supergravity (SUGRA) mediated SUSY 
breaking the soft SUSY breaking masses are generated at the Planck scale, and the low 
energy soft masses relevant for physical proceses such as LFV are therefore subject 
to radiative corrections in running from the Planck scale to the weak scale. Off- 
diagonal slepton masses in the SCKM basis can arise both directly at the high energy 
scale (due to the effective SUGRA theory which is responsible for them), or can be 
radiatively generated by renormalisation group running from the Planck scale to the 
weak scale, for example due to Higgs triplets in GUTs or right-handed neutrinos in 
see-saw models which have masses intermediate between these two scales. 

Neutrino experiments confirming the Large Mixing Angle (LMA) MSW solution 
to the solar neutrino problem j3] taken together with the atmospheric neutrino data 
[5] shows that neutrino masses are inevitable [B]- The presence of right-handed neu- 
trinos, as required by the see-saw mechanism for generating neutrino masses, will 
lead inevitably to LFV, due to running effects, even in minimal SUGRA (mSUGRA) 
which has no LFV at the high energy GUT or Planck scale Therefore, merely 

assuming SUSY and the see-saw mechanism, one expects LFV to be present. This 
has been studied, for example, in mSUGRA models with a natural neutrino mass 
hierarchy [5]. There is a large literature on the case of minimal LFV arising from 
mSUGRA and the see-saw mechanism 

Despite the fact that most realistic string models lead to a low energy effective 
non-minimal SUGRA theory, such theories have not been extensively studied in the 
literature, although a general analysis of flavour changing effects in the mass insertion 
approximation has recently been performed [TT]. Such effective non-minimal SUGRA 
models predict non-universality of the soft masses at the high energy scale, dependent 
on the structure of the Yukawa matrices. Moreover there can be additional sources 
of LFV which also enter the analysis. For example, realistic effective SUGRA theo- 
ries arising from string-inspired models will also typically involve some gauged family 
symmetry which can give an additional direct (as opposed to renormalisation group 
induced) source of LFV. This is because the D-term contribution to the scalar masses 
generated when the family symmetry spontaneously breaks adds different diagonal 
elements to each generation in the 'theory' basis, and generates non-zero off-diagonal 
elements in the SCKM basis, leading to LFV. This effect depends on the strength of 
the D-term contribution, which is expected to be close in size to the size of the un- 
corrected scalar masses. There can also be a significant contribution to non-universal 
trilinear soft masses leading to flavour violation ^21 EH E] arising from the F-terms 
of scalars associated with the Yukawa couplings (for example the flavons of Froggatt- 
Nielsen theories). 

The purpose of the present paper is to catalogue and quantitatively study the im- 
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portance of all the different sources of LFV present in a general non-minimal SUGRA 
framework, including the effects of gauged family symmetry. Although the different 
effects have all been identified in the literature, there has not so far been a coher- 
ent and quantitative dedicated study of LFV processes, beyond the mass insertion 
approximation, which includes all these effects within a single framework. In order 
to quantify the importance of the different effects it is necessary to investigate these 
disparate souces of LFV numerically, both in isolation and in association with one 
another, within some particular SUGRA model. To be concrete we shall study the 
effective SUGRA models of the kind considered in ^3] which have a sufficiently rich 
structure to enable all of the effects to be studied within a single framework. Within 
this class of models we shall consider specific benchmark points in order to illus- 
trate the different effects. Some of these benchmark points were already previously 
considered in [T^j. However in the previous study the important effect of D-terms 
arising from the Abelian family symmetry was not considered. Here we shall show 
that such D-terms are in fact calculable within the framework of the particular model 
considered, and can lead to significant enhancement (or suppression) of LFV rates, 
depending on the particular model considered. 

The outline of this paper is as follows. In Section|2]we discuss soft supersymmetry 
breaking masses in supergravity. In Section El we summarise the flavour problem and 
catalogue the distinct sources of lepton flavour violation. In Section 01 we outline 
the aspects of the specific models that we shall study. In Section we discuss the 
soft SUSY-breaking sector within this class of models, parameterise the F-terms, and 
write down the soft terms, including the D-term contribution to the scalar masses. 
We also discuss the D-terms associated with a family symmetry that are expected to 
lead to large lepton flavour violation. In Section |H1 we specify two models in this class, 
define benchmark points and present the results of our numerical analysis of LFV for 
these benchmark points. Section [7| concludes the paper. 



2 Soft terms from supergravity 

We summarise here the standard way of getting soft SUSY breaking terms from 
supergravity. Supergravity is defined in terms of a Kahler function, G, of chiral 
superfields (<fr = h,C a ). Taking the view that supergravity arises as the low energy 
effective field theory limit of a string theory, the hidden sector fields h are taken to 
correspond to closed string moduli states (h = S,Ti), and the matter states C a are 
taken to correspond to open string states. In string theory, the ends of the open 
string states are believed to be constrained to lie on extended solitonic objects called 
Dp— branes. 

Using natural units, 
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K(4>, 4>) is the Kahler potential, a real function of chiral superflelds. This may be 
expaned in powers of C a : 



K — K(h, h)+Ka b (h, h)CaC b 



-Z ab (h, h)C a C b + h.c. 



+ ... 



(2) 



Ka b is the Kahler metric. W(4>) is the superpotential, a holomorphic function of 
chiral superflelds. 

W = w(h) + lfi ab (h)C a C b + l -Y abc C a C b C c + ... (3) 

We expect the supersymmetry to be broken; if it is broken, then the auxilliary 
fields ^ for some <fi. As we lack a model of SUSY breaking, we introduce 
Goldstino angles as parameters that will enable us to explore different methods of 
breaking supersymmetry. We introduce a matrix, P that canonically normalises the 
Kahler metric, P^KjjP = 1 2 ^3]. We also introduce a column vector which 
satisfies 0^0 = 1. We are completely free to parameterise in any way which 
satisfies this constraint. 

Then the un-normalised soft terms and trilinears appear in the soft SUGRA break- 
ing potential [TH] 

Vsoit = ml b CaC b + (^A abc Y abc C a C b C c + h.c^j + ... (4) 

The non-canonically normalised soft trilinears are then 
W 



\W\ 



K m Y abc + d m Yabc- ((k ^ OrnKeaYdbc 

+ (a^b) + (a^ c) )1 (5) 



In this equation, it should be noted that the index m runs over h, C. However, 
by definition, the hidden sector part of the Kahler potential and the Kahler metrics 
is independent of the matter fields. 

Assuming that the terms dcY abc ^ 0, the canonically normalised equation for the 
trilinear is 

A abc = Fj [if r - dj In (k a k b K^ ] + F m d m In Y abc (6) 

If the Yukawa hierarchy is taken to be generated by a Froggatt-Nielsen (FN) field, 
0, such that Y oc P , then we expect F^ oc m 3 / 2 (f), and then F^c^lnF oc m 3 / 2 , and 
so even though these fields are expected to have heavily sub-dominant F-terms, they 
contribute to the trilinears on an equal footing as the moduli. 

2 The subscripts on the Kahler potential Kj means di K. However, the subscripts on the F-terms 
are just labels. 
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If the Kahler metric is diagonal and non-canonical, then the canonically normalised 
scalar mass-squareds are given by 

ml^ml/.-FjFjdjdj (in K a ), (7) 
and the gaugino masses are given by 

where f a is the 'gauge kinetic function', a enumerates D-branes in the model. In 
type I string models without twisted moduli these have the form f 9 = S, / 5i = Tj. 

Specifically, we use a Kahler potential that doesn't have any twisted-moduli |17j : 

K = -\n(s + S -\C^\ 2 -\Cl 2 \ 2 ^) -In (Vi +T X - jCt) 2 - jCfl 

-In (t 2 -T 2 - \Ci\ 2 - |C 3 5l | 2 ) - In (t 3 - T 3 - |Cf| 2 - |C 2 51 | 2 - |Cf \ 2 

|2 



| C 5 l52 |2 ^ | C95l 



(5 + S) 1/2 (T 3 + f s ) 1/2 (T 2 + T 2 ) 1/2 (T 3 + T 3 ) 1/2 

| C 95 2 |2 



(Ti + Ti) 1 ^ 2 (T 3 + T 3 ) 1/2 



(9) 



The notation is that the field theory scalars, the dilaton S and the untwisted 
moduli Ti originate from closed strings. Open string states C\ are required to have 
their ends localised onto D— branes. The upper index then specifies which brane(s) 
their ends are located on, and if both ends are on the same brane, the lower index 
specifies which pair of compacitified extra dimensions the string is free to vibrate in. 



3 Sources of lepton flavour violation 

There are two parts to the flavour problem. The first is understanding the origin of the 
Yukawa couplings (and heavy Majorana masses for the see-saw mechanism), which 
lead to low energy quark and lepton mixing angles. In low energy SUSY, we also need 
to understand why flavour changing and/or CP violating processes induced by SUSY 
loops are so small. A theory of flavour must address both problems simultaneously. 
For a full discussion of this see the review |3] . 

There are two contributions that can lead to large amounts of flavour violation. 
The first is the non-alignment of the trilinear soft coupling matrices to the corre- 
sponding Yukawa matrices, due to the contribution F m d m Y, m = {H, H, 6, 8}. The 
reasons why this can lead to large flavour violation are have been given before |14j . 
where a numerical investigation of a model very similar to those considered herein 
finds that there is a large amount of flavour violation. The second contribution can 



4 



come from scalar mass matrices which are not proportional to the identity in the 
theory basis, and lead to off-diagonal entries in the SCKM basis, resulting in flavour 
violation. 

In this section we begin by defining the SCKM basis, and in the following sub- 
sections we systematically discuss a number of distinct sources of Lepton Flavour 
Violation (LFV) in SUGRA models. As well as considering generic SUGRA models, 
we also allow for a family symmetry, which easily lead to non-universal scalar mass 
matrices, and non-aligned trilinear matrices 3 . 



3.1 The SCKM Basis 

The most convenient basis to work in for considering flavour violating decays, such 
as /i — > ey is the super-CKM (SCKM) basis, which is the basis where the Yukawa 
matrices are diagonal. If we define the unitary rotation matrices Uf,Vf by 

yl = u f yfv h ( io ) 

such that has positive eigenvalues. To convert the physical mass matrices to the 
SCKM basis, we rotate by the relevant matrix; Uf for the left-handed scalar matrices 
and Vf for the right-handed. Then flavour violation is proportional to the off-diagonal 
elements in the SCKM basis, and is suppressed by the diagonal values. The selectron 
mass matrix is 6 by 6, and the sneutrino mass matrix is 3 by 3 4 . The selectron mass 
matrix is 

/ Y e Y e h>l + W(gi-gj)l + m 2 LL -Y e v 2 /j + A e v x 

Tft — I + 

E \ -Y et v 2 fi* + A e vx Y^Y e v\ + \v 2 g\\ + m\ 

where v 2 = v\ — vf. The sneutrino mass matrix is then 

ml={Y^vl + \v 2 gll + ml^ (12) 

Off diagonal elements in any of the 3 by 3 submatrices in the SCKM basis will 
lead to flavour violation. We will now consider the LL block of m|,. The arguments 
follow for any other block of m| or m~. The transformation to the SCKM basis is 
carried out by 

m\ LL = U e (y*Y*v{ + \v\g\ - g\)\ + Xj\ (13) 

U e is unitary, and U e Y e V^ is diagonal, so the first two terms will be diagonal. Any 
off-diagonality must come from the third term. If this is proportional to the identity at 

3 By non-aligned trilincars, wc mean that A^jY^ ^ constant (no sum). 

4 The heavy right handed neutrinos cause the right-handed part of the sneutrino mass matrix to 
decouple by the electroweak scale. 
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the GUT scale, it will be approximately equal to the identity at the electroweak scale, 
which is the scale we should be working at. The fact that this is only approximate 
is due to the presence of the right handed neutrino fields in the running of the soft 
scalar mass squared matrices. If, however, the soft mass squared matrices are not 
proportional to the identity at the GUT scale, then large off-diagonal values will 
be generated when rotating to the SCKM basis, unless the rotation happens to be 
small. Generally this won't be the case. Since the family D-term contribution is not 
proportional to the identity 5 this will usually be the case 6 and so we expect large 
flavour violation in models with Abelian family symmetries when the D-terms correct 
the scalar mass matrices. 

3.2 The relevance of the Yukawa textures 

There is one subtlety concerning the size of the off-diagonal elements of the scalar mass 
matrices in the SCKM basis. This comes back to the definition of the SCKM basis 
as the basis in which the Yukawa matrices are diagonal. The larger the SCKM trans- 
formation between any 'theory' basis and the mass eigenstate basis for the Yukawa 
matrices, the larger the SCKM transformation that must be performed on the scalar 
mass matrices in going to the SCKM basis, hence the larger the off-diagonal elements 
of the scalar mass matrices in the SCKM basis generated from non-equal diagonal 
elements in the 'theory' basis. 7 The larger the off-diagonal entries in the SCKM basis 
compared to the diagonal ones, the greater will be the flavour violation. Also, the 
greater the mass difference between the diagonal elements in the 'theory' basis, the 
greater the size of the off-diagonal entries produced when rotating from the 'theory' 
basis, hence the larger the flavour violating effect. Clearly these effects are sensitive 
to the size of the transformation required to go to the SCKM basis, which in turn 
is sensitive to the particular choice of Yukawa textures in the 'theory' basis. In this 
way, the choice of Yukawa texture can play an important part on controlling the 
magnitude of flavour violation, and we shall see examples of this later. 

3.3 Running effects 

Consider the case where, at the high-energy scale, the scalar mass matrices are pro- 
portional to the identity matrix and each soft trilinear coupling matrix is aligned to 

5 This statement assumes that the generational charges are not the same for both left- and right- 
handed fields. This would remove the point of the family symmetry generating the fermion mass 
hierarchy. 

6 One can, however, imagine some model with aberrant points in its parameter space where a 
non-universal non-zero D-term corrects a non-universal base mass matrix to give a universal net 
mass matrix. 

7 Note that the D-terms make us sensitive to right-handed mixings in the Yukawa matrices, so the 
non-universal family charge structure for the right-handed scalar masses may lead to a non-universal 
generational hierarchy in the right-handed scalar mass matrices. 
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the corresponding Yukawa matrix: 

("'?), ; <r\! , (A^^AlYl (14) 

This is often referred to as mSUGRA. In the quark sector, due to the quark flavour 
violation responsible for CKM mixing, when the scalar squark mass matrices are run 
down to the electroweak scale, they will run to non-universal scalar mass matrices and 
non-aligned trilinear coupling matrices. If this is the case, then in the SCKM basis, 
which is the basis where the Yukawa matrices are diagonal, off-diagonal elements in 
the scalar squark mass matrices or the trilinear squark mass matrices lead to flavour 
violation. 

In the lepton sector, in the absence of neutrino masses the separate lepton flavour 
numbers are conserved and mSUGRA will not lead to any LFV induced by running the 
matrices down to low energy. However, in the presence of neutrino masses, with right- 
handed neutrino fields included to allow a see-saw explanation of neutrino masses and 
mixing angles, the separate lepton flavour numbers will be violated and, even in the 
mSUGRA type scenario, running effects will generate off-diagonal elements in the 
scalar mass matrices in the SCKM basis, resulting in low energy LFV. 

3.4 Diagonal scalar mass matrices not propotional to the unit 
matrix 

3.4.1 Non-minimal SUGRA 

In non-minimal SUGRA the scalar mass matrices may be diagonal at the high-energy 
scale, but not proportional to the identity. In this case, there will be non-zero off- 
diagonal elements in the SCKM basis even with no contribution to running effects, 
or contribution from the trilinear coupling matrices. 

One way of getting diagonal mass matrices not proportional to the unit matrix 
is from a SUGRA model corresponding to the low energy limit of a string model 
with D-branes. If each generation from the field theory viewpoint corresponds to a 
string attaching to different branes, then the masses predicted in the SUGRA can be 
different. This leads to diagonal but non-universal scalar mass matrices. 

3.4.2 D-term contributions from broken family gauge groups 

Another way of getting diagonal mass matrices not proportional to the unit matrix 
is by having a model with a gauge family symmetry, which is broken spontaneously. 
When the Higgs which breaks the family group, the flavon, gets a vev, it contributes 
a squark (slepton) mass contribution through the four point scalar gauge interaction 
which has two flavons and two squarks (sleptons). 

To make the point more explicitly, consider a U(l) family group. Then the mass 
contribution is proportional to the charge under the family symmetry. As the point of 
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a family symmetry is to explain the hierarchy of fermion masses, small quark mixing 
angles and large neutrino mixing angles, the charges are usually different. 

Then, even if the mass matrix starts off as a universal matrix, it will be driven 
non-universal by the D-term contribution: 



Din 



<?L1 



ffl, 



mf 



D 



QL2 



(15) 



11)7, 



<?L3 



3.5 Non-aligned trilinears 
3.5.1 Non-minimal SUGRA 

One way of getting non-aligned trilinear matrices is by having the same sort of non- 
minimal SUGRA setup that leads to diagonal but non-universal mass matrices, as 
described in Section 13.4.11 From the supegravity equations from Section the tri- 
linears that appear in the soft Lagrangian, Ay will be non-aligned if the trilinears 
predicted by the SUGRA model, Ay are not democratic, i.e. if Ay ^ constant. From 
a string-inspired/SUGRA standpoint, if each generation is assigned to a different 
brane and extra-dimensional vibrational direction, then in general we expect A a t, c to 
be different, due to the differing values of the Kahler metrics K a for the different brane 
assignents Cj. When is transformed to the SCKM basis at the electroweak scale, 
there will then be large off-diagonal elements which contribute to flavour violating 
processes. 



3.5.2 Flavon contributions from the Yukawa couplings 

In general when one considers a family symmetry in order to understand the origin 
of the Yukawa couplings, the new fields arising from this can develop F-term vevs, 
and contribute to the supersymmetry breaking F-terms in a non-universal way. This 
leads to a dangerous source of flavour violating non-aligned trilinears: [T2 | IT3 t lT4 ^ ll8j. 

AA = Fgdg In Y (16) 

where the Yukawa coupling Y in Eq. (JTfij) arises from the an effective FN operator 
and is a polynomial of the FN field 9, Y ~ 6 n , leading to 

71 

AA = F e d e ln6 n = F e -. (17) 

a 

However the auxiliary field is proportional to the scalar component, 

F e oc m 3 / 2 9 =>• AA oc nm 3 / 2 . (18) 

An example of this with an arbritary U(l) family symmetry is 

/ q y(i,j) 

Yij = aij ( j ~ rn 3/2 p(z,j) (19) 
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5i-brane 




Figure 1: A generic type I string construction involving two sets of perpendicular D5-branes 
embedded within a D9-brane, where the D5-brane world-volumes intersect at the origin. 
Charged chiral fields appear as open strings with both ends attached to the same D-brane 
Cj 1 and Cj, or different branes C 5lS2 and C 95 *. Closed strings (S,Ti) can live in the full 
lOd space, although orbifolding leads to closed strings (twisted moduli Y^j localised at 4d 
fixed points within the Z?5j-brane world- volume. 

The arbiratry couplings, all of which should be 0(1) for the symmetry to 

be considered natural. The p(i,j) are integers appearing as a power for the ij-th 
element of the above Yukawa, and it comprises the sum of the family charges for the 
ith-generation left-handed field and jth-generation right-handed field. In principle, 
if the Yukawa texture is set up so that each power is different, then each element 
in Aij will be different from each other, and the physical trilinear matrix, Ay will 
be non-aligned to the corresponding Yukawa. Due to the dependence on the charges 
of the different fields, this contribution to the trilinears is not diagonalised when we 
transform to the SCKM basis. 

4 Intersecting D-brane models with an Abelian 
family symmetry 

4.1 Symmetries and symmetry breaking 

In order to study the effects of LFV elucidated in the previous section, it is necessary 
to specialize to a particular effective non-minimal SUGRA model which addresses 
the question of flavour (i.e. provides a theory of the Yukawa couplings). The specific 
model we shall discuss is defined in Table ^ This model is an extension of the 
Supersymmetric Pati-Salam model discussed in ref.^H], based on two D5 branes 
which intersect at 90 degrees and preserve SUSY down to the TeV energy scale. The 
generic D-brane set-up that we use is illustrated in Figd where the string assignment 
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notation is defined. The gauge group of the 5i sector is 11(4)^ x U(2)^ x 1/(2)^, and 
the gauge group of the 52 sector is 11(4)^ (e.g., we assume the U(2)l^r of the 52 sector 
are broken). The symmetry breaking pattern of this model takes place in two stages, 
which we assume occur at very similar scales ~ Mx- In the first stage, the U(4) 
groups are broken to the diagonal subgroup via diagonal VEV's of bifundamentals; 
the resulting theory is an effective Pati-Salam model (with additional [/(l)'s) which 
then breaks to the MSSM (and a number of additional U(l)'s) via the usual Higgs 
pair of bifundamentals. The string scale is taken to be equal to the GUT scale, about 
3 x 10 16 GeV. 

The symmetry breaking pattern leads to the following relations among the gauge 
couplings of the SM gauge groups in terms of the gauge couplings g 5l and g 52 associ- 
ated with the gauge groups of the 5i and 5 2 sectors: 

95i 952 



9s 
9i 



gl + 9l 2 



95 1 = 92R 

>/%302 



9a 



2g 2 2 



(20) 

(21) 
(22) 



Field 


U(4)W 


U(2)g> 


U(2)£ j 


U(4)( 2 ) 


Ends 


U(l)j? charge 


h 


1 


2 


2 


1 









4 


2 


1 


1 




QL3 


F 3 


4 


1 


2 


1 




QR3 


F 2 


1 


2 


1 


4 


(^5i5 2 


QL2 


F 2 


1 


1 


2 


4 




QR2 


F X 


1 


2 


1 


4 


^5i5 2 


Qli 


F X 


1 


1 


2 


4 


^5i5 2 


Qri 


H 


4 


1 


2 


1 


cf 


Qh 


H 


4 


1 


2 


1 


c 2 51 


-qn 




4 


1 


1 


4 


(J5i52 




¥2 


4 


1 


1 


4 


(^5i5 2 






6 


1 


1 


1 








6 


1 


1 


1 






9 


1 


1 


1 


1 




-1 


9 


1 


1 


1 


1 


All 


1 



Table 1: The particle content of the 42241 model, and the brane assignments of the 
corresponding string. Note that the string assignment of 6 is allowed to be any of 
{C 5l5;2 , C^ 1 , Cf 1 , Cf 1 }, giving a slightly different model in each case. The U(1)f charges 
will change, also generating a slightly different model in each case. We will return to this 
in more detail in Sectional 

The extension is to include an additional U(1)f family symmetry and the FN 
operators responsible for the Yukawa couplings as in [20] (see also [2Ij)- The charges 
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under the Abelian symmetry U(1)f are left arbritary for now. The present '42241' 
Model is the same as the model considered in J3], with the following modifications 
considered; firstly, we allow the Froggatt-Nielsen field 9 to be either an intersection 
state or attached to the 5i brane. The location of 9 dramatically changes the value 
of the D-term contribution to the scalar masses coming from the FN sector. 

The quark and lepton fields are contained in the representations F, F which are 
assigned charges under U(l)jr. In Table [Owe list the charges, string assignments and 
representations under the string gauge group U(A)^ ® U(2)^ ® U(2)^ ® U{A) <y2 \ 

The field h represents both Electroweak Higgs doublets that we are familiar with 
from the MSSM. The fields H and H are the Pati-Salam Higgs scalars; 8 the bar on 
the second is used to note that it is in the conjugate representation compared to the 
unbarred field. 

The extra Abelian U(l)p gauge group is a family symmetry, and is broken at the 
high energy scale by the vevs of the FN fields [22] 9,9, which have charges —1 and +1 
respectively under U(l)p- We assume that the singlet field 9 arises as an intersection 
state between the two D5— branes, transforming under the remnant U(l)s in the 4224 
gauge structure. In general the FN fields are expected to have non-zero F-term vevs. 

The two SU (4) gauge groups are broken to their diagonal subgroup at a high scale 
due to the assumed vevs of the bifundamental Higgs fields ipi, ip 2 |T5] . The symmetry 
breaking at the scale M x 

SU(4) <g> SU(2) L <g> SU(2) R -> SU(3) ® SU(2) L ® U(l)y (23) 

is achieved by the heavy Higgs fields H, H which are assumed to gain vevs [21] 

(H ab ) = (u H ) = V525 b 2 ~ M x ; (H ax ) = (V H ) = V5 4 J 2 X ~M X (24) 

This symmetry breaking splits the Higgs field h into two Higgs doublets, hi, h 2 . Their 
neutral components then gain weak-scale vevs. 

(h° 1 )=v 1 ; (h° 2 )=v 2 ; tfw/3 = v 2 /v 1 . (25) 

The low energy limit of this model contains the MSSM with right-handed neutrinos. 
We will return to the right handed neutrinos when we consider operators including 
the heavy Higgs fields H , H which lead to effective Yukawa contributions and effective 
Majorana mass matrices when the heavy Higgs fields gain vevs. 



4.2 Yukawa operators 

The (effective) Yukawa couplings are generated by operators involving the FN field 9 
with the following structure: 9 |21| : 



fHH\ n ( 9 



P(i,j) 



8 We will also refer to these as "Heavy Higgs"; this has nothing to do with the MSSM heavy 
neutral higgs state H°. 

9 The field 6 will not enter the Yukawa operators because FiFjh will be positive for any i,j. 
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where the integer p(i, j) is the total U(1)f charge of Fi + Fj + h and HH has & U(1)f 
charge of zero. The tensor structure of the operators in Eq. f)26|) is 

{0)^° z = F aa Fp x hlH lz H? w e^ (27) 

One constructs SU(4)ps invariant tensors that combine 4 and 4 representations 
of SU(4)p5 into 1, 6, 10, 10 and 15 representations [21] • Similarly we construct 
SU(2) R tensors R*^ that combine the 2 representations of SU(2) into singlet and 
triplet representations. These tensors are contracted together and into 0°^™ z to 
create singlets of SU(4)ps, SU(2) L and SU(2)p. Depending on which operators are 
used, different Clebsch-Gordan coefficients (CGCs) will emerge. 

We will return to these in section El when we define the two models that we will 
be using for the numerical analysis. 



4.3 Majorana operators 

We are interested in Majorana fermions because they can contribute neutrino masses 
of the correct order of magnitude via the see-saw effect. The operators for Majorana 
fermions are of the form 



There do not exist renormalisable elements of this infinite series of operators, so 
n < 1 Majorana operators are not defined, except in the (3, 3) element. We assume 
that a (3, 3) neturino mass term is allowed at leading (but non-renormalisable) order. 
A similar analysis goes through as for the Dirac fermions; however the structures only 
ever give masses to the neutrinos, not to the electrons or to the quarks. 



5 Soft supersymmetry breaking masses 
5.1 Supersymmetry breaking F-terms 

In [HI] it was assumed that the Yukawas were field-independent, and hence the only 
F-vevs of importance were that of the dilaton (S), and the untwisted moduli {T l ). 
Here we set out the parameterisation for the F-term vevs, including the contributions 
from the FN field 8 and the heavy Higgs fields H, H. Note that the field dependent 
part follows from the assumption that the family symmetry field, 9 is an intersection 
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state. 

F s = V3m 3/2 (S + S)X s (29) 
F Ti = V3m 3/2 (Ti + Ti) X Ti (30) 

Fjjab = VZm z/2 H ab (S + Sy X H (31) 

F-jj ax = V3m 3/2 H ax {T 3 + T 3 yX w (32) 

F e = V3m 3/2 e(S + Sy(T 3 + T 3 ^X e (33) 

We introduce a shorthand notation: 

F H H = F H « b H ab ; FffH = ^ F^JS^. (34) 

ab ax 

The F-terms above use values of S and Tj which are given in terms of the gauge 
couplings as: 

Air 4-7T 
Re(S) = — Re(Ti) = — . (35) 

The gauge couplings g^g^ are given from Eq. (|20 J) -(|23 J) [TH] as 

9bi = 92, 95 2 = i 9 9 , (36) 

V #2 - #3 

where we shall assume that at the scale Mx we have, 

g 2 = 0.7345, g 3 = 0.6730. (37) 
The values of g$, g$ 3 are assumed to be equal and are obtained from the string relation 

327T 2 C^^j = 0906i06a05 3 > (38) 

as 

(7 9 = g, 3 = 0.0266, (39) 



where we have taken 

fiM =2.77xl0~ 6 . (40) 
\M Pl ) 

These rather small gauge couplings imply 

Re(5) = Re(T 3 ) = 0.877. (41) 

In [T3| the string relation was not used and it was assumed incorrectly that g 9 = 
g$-i — g2 which resulted in Re(5) = 27.7. 
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5.2 Soft scalar masses 



There are two contributions to scalar mass squared matrices, coming from SUGRA 
and from D-terms. In this subsection we calculate the SUGRA predictions for the 
matrices at the GUT scale, and in the next subsection we add on the D-term contri- 
butions. 

The SUGRA contributions to soft masses are detailed in Section [21 
From Eq. (JJJ) we can get the family independent form for all scalars: 



where 



m, 



Til 



3/2 



(42) 







a 






m R = 


m 3/2 


a 




(43) 






h R . 






m h = 


m 3/ 2 ( 


1 - 3X1) 




(44) 


m H = 


ml /2 (l-3X 2 s ) 




(45) 


m TT = 


m l/2i 


1 - 3X 2 Ts ) 




(46) 


m 2 e = 


m l/2 






(47) 


a 


= 1 - 


-\{xl + x^) 


(48) 


b L 


= 1 - 


- 3X^ z 




(49) 


b R 


= 1 - 


- 3Xf. 2 




(50) 



Qi 



and mh r . m 



Here m\ represents the left handed scalar mass squared matricies m 
represents the right handed scalar mass squared matricies fnfj R , vt? Dr ^ m \ R and ™ 2 n r - 
A discussion of the equations for m| can be found in Section 15.41 



5.3 D-term contributions 

There are two D-term contributions to the scalar masses. The first is the well known 
[2*3*1 E] contribution from the breaking of the Pati-Salam group to the MSSM group. 
Note that these D-terms are different to those quoted in the references above as we 
now consider the D-terms generated by breaking a family symmetry. See Appendix iBl 
for a full derivation. The second D-term comes solely from the breaking of the U(l) 
family symmetry. The corrections lead to the following mass matrices: 
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qLi 



m 



Qi 



m. 



m 



in 



m, 



m 



N R 



m, 



m h 



m 2 L + l(gi)D-j I + | q L2 \g 2 F D 2 
m 2 L -l(3g 2 )D 2 H + \ q L2 I g 2 F D 2 



m 2 R -l(gi-2g 2 2R )D 2 H 



qm 



| 9 2 F D 2 

qm, 

q m g 2 F D 2 , 

qm) 

qm \ 

m\ + l{3g 2 - 2g 2 m )D 2 H + q m g 2 F D 2 

qm) 

qm g 2 F D 2 e 

qm 



qm 



m R +l(3gi + 2g 2 2R )D 2 H 



ml-2g 2 2R D 2 H , 
m 2 hl + 2g 2 2R D 2 H , 



(51) 
(52) 
(53) 
(54) 
(55) 

(56) 

(57) 
(58) 



The charges quARj are the charges under U{l)p of Fj and Fj respectively, as 
shown in Tabled The correction factors D 2 ,, D 2 H are calculated explicitly in Appendix 
IHlin terms of the gauge couplings and soft masses as 10 

D h = Ag 2 R \Qg 2 ^ ~ + qH ^ e ~ ^ 



1 2 

mi — m- r 



D e = tt~o (60) 



2g 



We note that the factors of g 2 F appearing in the mass matrices are cancelled by 
the 4- in the definition of Do. 

9f 

The D-terms associated with the family symmetry depend on the charges of the 
left-handed and right-handed matter representations F, F under the family symme- 
try. It is well known 11 that for Pati-Salam, one can choose any set of charges, and 
there will be an equivalent, shifted set of charges that are anomaly free due to the 
Green-Schwartz anomaly cancellation mechanism. The charges used for the D-term 
calculation should be the anomaly free charges. 

10 qn is defined to be — qjj 3 , thus qn = | for Model 1 and qn = 1 for Model 2. See Appendix iBl 
for more details. 

11 For an explanation, see for example 24 . 
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The gauge couplings and mass parameters in Eqs. (|5*51 , 15U|) are predicted from 
the model, in terms of the X parameters and m 3 / 2 as shown in Eqs. ([451 - 117)1 and 
Eqs. (J52l - Note that the D-terms will be zero if Xs = Xt v or if the 9 brane 
assignment is the same as 9. Choosing the second of these conditions is useful since 
it gives a comparison case where there are no U(l)p D-terms; this comparison will 
make the D-term contribution to flavour violation immediately apparent. 

5.4 Magnitude of Z^-terms for different 9 assignments 

The main point worth emphasising is that in the string model the magnitudes of 
the D-terms are calculable. We have assumed throughout that the FN field 9 is an 
intersection string state C 5l5i , but have not specified the string assignment of 9. Thus 
m| takes various values depending on the string assignment for 9. 

From Eq. (j6T?j) . we see that we have calculable D-terms, 

2g 2 F D 2 e = - m§ , (61) 

so the value of D 2 e depends on the choice of where the 9 field lives. We use Table Q 
and Eqs. (j4"21 - loTH) to quantify the LVterm for each possible 9 string assignment. As 
9 always lives at the intersection, on C 5l5i , our first choice of 9 on C 5lb2 is trivial: it 
gives D 2 = 0. For 9 on Cf 1 , m| is equivalent to m|, as this is also on C^ 1 . So using 
Eq. ()44|) for m| and Eq. (|47J) for m 2 B in Eq. (|61|). we have 

Cf : 2g 2 F D 2 d = ^m 2 3/2 (X 2 -X 2 3 ). (62) 
Similarly, the other two choices yield 

: 2g 2 F D 2 = ~m 2 /2 (X 2 - X 2 3 ) (63) 
C3 51 : 2g 2 F Dl = - 3 -m 2 3/2 (X 2 -X 2 2 ). (64) 

In this paper, all models use X F3 = X F2 = X Fi = X T , so Eqs. (p)3*|) and (j6^j) are 
equal to each other, and opposite in sign to Eq. (jH2|)- 

5.5 Soft gaugino masses 

The soft gaugino masses are the same as in JH] , which we quote here for completeness. 
The results follow from Eq. flSJ) applied to the SU(4) <8> SU(2) L <g> SU(2) fl gauginos, 
which then mix into the SU(3) ® SU(2)l (8) U(l)y gauginos whose masses are given 
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by 



M 3 
M 2 
Mi 



m 



3/2 



(Ti+Ti) + (T 2 + T 2 ) 
y/3m 3/2 X Tl 

V3m 3/2 

Ut 1 + t 1 ) + Ut 2 + t 2 ) 



=- [(Ti + TOX Tl + (T 2 + T 2 )X T2 ] 



!(T 1 + T 1 )X Tl + ^(T 2 + T 2 )X T2 



(65) 
(66) 
(67) 



The values of T\ + Ti and T 2 + T 2 are proportional to the brane gauge couplings g 5l 
and g§ 2 , which are related in a simple way to the MSSM couplings at the unification 
scale. This is discussed in [T9*] . 

When we run the MSSM gauge couplings up and solve for g 5l and g 52 we find that 
approximate gauge coupling unification is achieved by T\ + T\ 3> T 2 + T 2 . Then we 
find the simple approximate result 



M 1 ^M 3 ^M 2 = V3m 3/2 X Tl . 



(68) 



5.6 Soft trilinear couplings 

So far we have considered the soft masses and the gaugino masses. The gaugino masses 
are the same as in ^H]- The soft masses have had both D-term contributions added 
onto the base values from The contributions to the soft masses and gaugino 

masses from the FN and heavy Higgs auxiliary fields is completely negligible due to the 
small size of their F-terms. However for the soft trilinear masses these contributions 
are of order 0(m 3 / 2 ) despite having small F-terms, so FN and Higgs contributions 
will give very important additional contributions beyond those considered in |19j . 

From Section |21 we see that the canonically normalised equation for the trilinear 

is 

A abc = Fj pCj - dj In {k a K b k)j ] + F m d m In Y abc (69) 

This general form for the trilinear accounts for contributions from non-moduli F- 
terms. These contributions are in general expected to be of the same magnitude as the 
moduli contributions despite the fact that the non-moduli F-terms are much smaller 
[25] . Specifically, if the Yukawa hierarchy is taken to be generated by a FN field, 9 such 
that Yij ~ 6 PiJ , then we expect Fg ~ tti 3 /2&, and then AAij = FgdglnYy ~ Vij m zi2 
and so even though these fields are expected to have heavily sub-dominant F-terms 12 
they contribute to the trilinears at the same order 0(m 3 / 2 ) as the moduli, but in a 
flavour off-diagonal way. 

Here we sum over m, which contains all the hidden sector fields: S, T, H, H,8. 



12 In our model the FN and heavy Higgs vevs are of order the unification scale, compared to the 
moduli vevs which are of order of the Planck scale. 
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In the specific D-brane model of interest here, the general results for soft trilinear 
masses, including the contributions for general effective Yukawa couplings are given 
in Appendix El From Eqs. (|26 [ ET j) we can read off the effective Yukawa couplings, 

Y hFF hFF ee ^)%{r)%H lz H<™e* h\F™F %. (70) 



Y —P x 
1 hFFuy 

Note the extra group indices that the effective Yukawa coupling Y hFF ^ has, and 
proper care must be taken of the tensor structure when deriving trilinears from a 
given operator. We can write down the trilinear soft masses, A, by substituting the 
operators in Eq. (jC.lJl into the results in Appendix |XJ Having done this we find the 
result: 13 



A = V3m 3 / 2 
where 



d\ + d H + p(i, j)de d\ + d H + p(i, j)dg d 2 + d H + p(i, j)d e 
di + d H + p(i,j)d e d 1 + d H + p(i,j)d e d 2 + d H + p(i, j)d e 
d 3 + d H +p(i,j)d e d 3 + d H +p(i,j)d e d 4 



(71) 



(72) 
(73) 

(74) 

(75) 
(76) 
(77) 

These results are independent of which string assignment we give to 9, since this 
field does not enter the Yukawa operators. 



d\ = 


X s — Xtx — X T2 


d 2 = 


—Xq — At 1 . — — Att, 

2 1 2 2 
-X s - X Tl - X T , 2 + -A Ts 


d 3 = 


d^ = 


— Atj 


d H = 


(S + S^XH + iTs + T^Xx 


d 9 = 


(S + S) l *(T 3 + T 3 )*X e 



6 Results 
6.1 Two models 

We will study two different models which have the same Yukawa textures, but different 
U(1)f charge structures. The anomaly-free [23 family charges are laid out in Tabled 
This first model is essentially the model studied in [201 E] , but with an extra operator 
in the (1,2) and (1,3) Yukawa matrix element to allow a non-zero Yf 2 and Yf 3 . The 
second model is defined such that all of the charges of left-handed matter are the same, 
causing the U(l)p D-term to left handed scalar mass matrices to not lead to extra 

13 We assume that a (3,3) Yukawa coupling appears at renormalisable order. This is why the ^33 
doesn't include contributions from du and dg. 
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flavour violation. This choice is a small change to the model, since two of the 'left- 
handed' charges are the same anyway, and is made because normally the left-handed 
contribution dominates over the right-handed contribution. This can be achieved 
by changing the order of the operators in the Yukawa textures and the arbritary 
couplings a, a', ... to compensate the change in the charge structure, as discussed in 
Appendix O Most of the results presented will be for Model 1. 

Model 2 differs from Model 1 in the charges under the Abelian family symmetry, 
and the compensating changes to ensure the same Yukawa textues. As discussed in 
Appendix O there are two expansion parameters in these models which we take to be 
equal e = 5 = 0.22. The powers of e in the first row of the Yukawa matrices are one 
lower, and we compensate that by increasing the powers of 5 in the first row. We do 
not have to shift the values of the a, a', a", a'" parameters since 5 = e. Were this not 
the case, we would have to shift the values by a factor of |. Model 2 is not meant to 
be a natural or realistic model, we use it as a tool to investigate the contribution to 
flavour violation from the U(l)p D-term correction to the left-handed scalar masses. 
The family charges are laid out in Table El 

The U(l)p charge of H must be equal and opposite to F 3 , and H must be the 
negative of this. This is due to the (3,3) element of the right-handed Majorana mass 
being allowed at leading order, so the U(l)p charges of F 3 and H must conspire to 
cancel for the operator of the Majorana fermions to be renormalisable. For Model 
1, the U(1)f charges of H and H are 5/6 and —5/6 respectively. For Model 2, the 
U(1)f charges of H and H are 1 and —1 respectively. One can use the relevent 
equations above to check that the anomaly coefficients do indeed satisfy the anomaly 
cancellation conditions. This gives us a different Dx-term for Model 2, as indicated 
by qn in Eq. (|59|) being different for Model 1 and Model 2. Note that D 2 B is the same 
in both models. 



Field 


Model 1 


Model 2 




Charge 


Charge 


Fx 


n 


1 


F 2 


8 

6 


1 


F 3 


5 
6 


1 


Fx 


19 
6 


3 


F 2 


7 
6 


1 


F 3 


5 
6 


-1 



Table 2: The family charges for Model 1 and Model 2. 

The values of the arbritary couplings are laid out in Appendix [U] in Table EJ This 
gives numerical values for the Yukawa elements which can be used in either model, 
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with the relevant values of Yf 2 and Yf 3 inserted instead of the texture zeros: 



Y U {M X ) 



Y d (M, 



Y e (M x ) 
Y U {M X ) 



2.159 x KT 06 

0.000 

0.000 

-1.661 x 10~ 04 
7.683 x 10~ 04 
-1.769 x 10~ 04 

-1.246 x 10~ 04 
1.537 x lO" 03 
-1.327 x 10~ 04 

2.159 x 10" 06 

0.000 

0.000 



5.606 x 10~ 04 
1.105 x 10" 03 
6.733 x 10~ 03 

-5.606 x 10" 04 
-5.343 x 10" 03 
3.133 x 10~ 02 

0.000 

2.432 x 10" 02 - 
3.133 x 10~ 02 

1.525 x 10~ 03 
8.290 x 10~ 04 
5.050 x 10~ 03 



5.090 x 10" 03 
0.000 

5.841 x 10~ 01 

1.018 x 10~ 02 
1.216 x 10~ 02 
3.933 x 10~ 01 

0.000 

-3.649 x 10~ 02 
5.469 x 10~ 01 

0.000 

3.923 x 10~ 01 
5.469 x 10~ 01 



(78) 
(79) 
(80) 
(81) 



The RH Majorana neutrino mass matrix for Models 1 and 2 has the numerical 
values: 

3.508 x 10 8 3.686 x 10 9 3.345 x 10 11 
3.686 x 10 9 8.313 x 10 10 5.886 x 10 12 (82) 
3.345 x 10 11 5.886 x 10 12 5.795 x 10 14 



M RR (M 



x, 



M- 



6.2 Benchmark points 

Since the parameter space for the models is reasonably expansive, and the intention 
is to compare different sources of LFV, it is convenient to consider five benchmark 
points, as follows. It should be noted that for all these points, we have taken all 
to be the same, = Xt, and also Xh = Xjj. X$ is taken to be zero throughout. 



Point 


X s 


Xt 


X H 


Xh 


X e 




A 


0.500 


0.500 


0.000 


0.000 


0.000 


0.000 


B 


0.535 


0.488 


0.000 


0.000 


0.000 


0.000 


C 


0.270 


0.270 


0.000 


0.000 


0.841 


0.000 


D 


0.270 


0.270 


0.595 


0.595 


0.000 


0.000 


E 


0.290 


0.264 


0.000 


0.000 


0.841 


0.000 



Table 3: Values of the X parameters for the five benchmark points, A-E. 



• Point A is referred to as "minimum flavour violation" . At the point Xs = X? 
the scalar mass matrices m 2 are proportional to the identity, and the trilinears 
A are aligned with the Yukawas. Also, if we look back to Eqs. (|59l60j) . (j43|) 
and (jH3), for Xs = Xt, which is the case for point A (and points C and D) 
we see that the value of both D-term contributions is zero. As such, both m 2 
and A would be diagonal in the SCKM basis in the absence of the RH neutrino 
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field, but in the presence of the see-saw mechanism off-diagonal elements are 
present in the SCKM basis leading to LFV. Since the four relevent soft masses 
are degenerate at this point, m3, 5l5o = m 2 5l = m 2 5l = m 2 5l , both D-term 

Cj/ Cj 1 Cg 1 

contributions are zero; D 2 e = = D\. 

• Point B is referred to as "SUGRA". With X$ 7^ Xt it represents typical flavour 
violation from the moduli fields, manifested as diagonal soft mass matrices not 
proportional to the unit matrix in the theory basis; this is the amount of flavour 
violation that would traditionally have been expected with no contribution from 
the Fh or Fq fields. This and point E are the only benchmark points investigated 
where Dg 7^ 0. 

• Point C is referred to as "FN flavour violation" . It represents flavour violation 
from the Froggatt-Nielsen sector by itself, arising as a non-alignment of the 
trilinear soft terms via the FN fields in the Yukawa operators, without any 
contribution to flavour violation from traditional SUGRA effects, since X$ = 
Xt as in point A. 

• Point D is referred to as "Heavy Higgs flavour violation" . It represents flavour 
violation from the heavy Higgs sector, arising as a non-alignment of the trilinear 
soft terms, via the Heavy Higgs fields in the Yukawa operators, without any 
contribution from either traditional SUGRA effects since Xg = Xt, or from FN 
fields since Fg = 0. 

• Point E combines features of points B and C, resulting in Froggatt-Nielsen 
flavour violation from X$, with SUGRA flavour violation from X$ 7^ Xt- This 
is the only point where we see the flavour violation from the Froggatt-Nielsen 
fields and the U(l) D-terms appearing at the same time. The numerical values 
for this point were obtained by taking the ratio from X$ and Xt for benchmark 
point B and applying it to benchmark point C. 

6.3 Varying Yf 2 and Yf 3 

Normally, the chargino contribution to LFV dominates. Since the Feynman diagram 
for this includes the left-handed sfermions, we would expect the D-term corrections 
to the left-handed slepton mass matrix to dominate the flavour violation. However, 
Model 2, as defined in Section \C2\ is set up to have universal left-handed charges, 
so the D-term correction from the breaking of U(1)f will not contribute to flavour 
violation (except that it will either add or remove some mass suppression). The D- 
term is limited in magnitude by the difference of m 2 , and m|, and although this is not 
a strong correction to the soft masses, it can contribute significantly to the lepton 
flavour violating branching ratios. 

The difference in \x — > e'j between Model 1 and Model 2 is negligible for Y± 2 = 
Yi 3 = 0. This should not be surprising, since the texture zero coming from Y( 2 = 
Y± 3 = will yield small mixing angles, resulting in small lepton flavour violation. 
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In order to get a picture of how great an effect the D-term contributions could 
have on the soft masses, it is necessary to examine a range of different values of Yf 2 
and Yi 3 . An extra operator contribution was added to the textures in Model 1 and 
Model 2, when compared to the model previously studied [Hj, to allow for variations 
of the order Y± 2 « 1CT 3 , for example. The gives 0(1) parameters a'[ 2 and a" 2 for 
Models 1 and 2 respectively. To be precise Yf 2 = 1.5 x 10~ 3 corresponds to a![ 2 and 
a'( 2 = 3.2, and Y* 3 = 1.5 x lCT 2 corresponds to a 13 = 0.320. 

6.4 Varying brane assignments for 9 

The 9 field is fixed to reside on the C 5l5i brane, but we allow the brane assignment 
of 9 to vary over the possibilities C\\ C\ x and Cf 1 . This gives us D-terms that 

are calculable in each case, rather than being free parameters. The assignment of 9 
to C 5l5i , which is the same assigmment as the 9 field, implies that the U(l)p D-term 
calculated in this case is zero. The other possibilities, given in Eqs. ()62l - 164|) . will 
highlight the contribution of the D-terms to lepton flavour violation. 

6.5 Numerical procedure 

The code used to generate all the data here was based on SOFTSUSY j2H], which 
is a program that accurately calculates the spectrum of superparticles in the MSSM. 
It solves the renormalisation group equations with theoretical constraints on soft 
supersymmetry breaking terms provided by the user. Successful radiative electroweak 
symmetry breaking is used as a boundary condition, as are weak-scale gauge coupling 
and fermion mass data (including one- loop finite MSSM corrections). The program 
can also calculate a measure of fine-tuning. The program structure has, in this case, 
been adapted to the extension of the MSSM considered in this paper. It is modified to 
include right-handed neutrino fields, and thus non-zero neutrino masses and mixing 
angles, generated via the SUSY see-saw mechanism. It is also set up to include 
the new D-term contributions considered herein, and to run over a series of string 
assignments for the 9 field. We use tan/3 = 50. 

Electroweak symmetry breaking provides a significant constraint on the results. 
The breakdown of electroweak symmetry breaking was responsible for the 'spike' 
feature that was shown in the plots for benchmark points A and B in [Hj. For 
the data above the spike, radiative electroweak symmetry breaking does not work 
properly, as the Z-boson mass becomes tachyonic. In the present paper such 'bad' 
regions where electroweak symmetry breaking fails are cut-off, however there is still 
a remnant of the spike left, which is why one can see a slight rise at the ends of the 
plots for our benchmark points A and B, as can be seen in Section 10)1 



22 



6.6 Numerical results 



We have now denned our two models, Model 1 and Model 2, and a set of five bench- 
mark points in Table El to examine within them. We have also set up what we will 
be varying apart from the gravitino mass in these models - the values of Y£ 2 , Y\z and 
the string assignment of 9 which gives different D-terms. We are now in a position to 
present our results. We shall focus on the branching ratios for \i — > c-f and r — > pq. 
The branching ratio for r — > c-f is not shown here as it does not constrain us beyond 
those limits placed by /i — > e^j and r — > pq. The experimental limit for r — > ej, at 
2.7 x 10~ 6 , is in fact far above the predicted rate for this process at all examined parts 
of the parameter space. In the following plots, we do not consider the 9 assignment 
to C3 1 as this is exactly the same as C2 1 , due to the degeneracy of the X Ti . Were we 
to allow the X^. to be non-degenerate, the phenomenological results of assigning 9 
to C2 1 and C3 1 would not be the same. The detailed spectrum will look different at 
each parameter point, but the general trend is for the physical masses to increase in 
magnitude as the gravitino mass increases. Thus too-high gravitino masses will start 
to reintroduce the fine-tuning problem resulting from the gluino mass being too high 
[2H], although we shall not discuss the detailed spectrum here. 

Figure 12 shows numerical results for BR(/i — > e^) for Model 1, plotted against the 
gravitino mass 7713/2, where each of the four panels (i) - (iv) correspond to each of the 
four benchmark points A - D. 14 Panel (i) of Figure 121 refers to benchmark point A, 
corresponding to minimum flavour violation, where the only source of LFV is from the 
see-saw mechanism, which for Model 1 is well below the experimental limit, shown as 
the horizontal dot-dash line. Panel (ii) of Figure El refers to benchmark point B, with 
LFV arising from SUGRA, with the FN and heavy Higgs sources of LFV switched off. 
In this case one can clearly distinguish the additional contributions to LFV arising 
from the D-terms. This makes benchmark point B the most phenomenologically 
interesting for the purposes of this study. The differing contributions stem from the 9 
string assignments, which are shown by the separate lines: C 5l5i (solid), Cf 1 (dashed), 
and C\ x (dot-dash). The C 5152 case shows the zero-D-term limit, where 9 and 9 are 
both intersection states, and hence conspire to cancel out Dg via their soft masses 
being degenerate, m? e — = 0. The other two locations for 9 then turn on the 
.tVterm contributions. With the D-terms switched on, Model 1 is experimentally 
ruled out over all parameter space shown here. Panel (iii) of Figure El refering to 
benchmark point C, is the Froggatt-Nielsen benchmark point, and for this case we see 
that the experimental limit is satisfied for 777.3/2 over 1400 GeV. Panel (iv) of Figure 121 
benchmark point D, shows the heavy Higgs point, for which the experimental limit 
is satisfied everywhere above 800 GeV. 15 

14 The solid line on each plot in Figure corresponds to the solid lines in Figure 1 of |14) . However 
there were errors in the code used to generate the previous data, and the corrected rates shown here 
differ to the previous results by up to two orders of magnitude. Furthermore, unlike the results in 
|14j . the results here do not exhibit a sharp spike for benchmark points A and B. For the data above 
the spike, radiative electroweak symmetry breaking does not work properly, as the Z-boson mass 
becomes tachyonic. This was not realized in the previous analysis. 

15 Note that the predictions in this figure are lower than the corresponding figure in |14j due to 
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Figure 2: Plots showing the Branching Ratio for /i — ► e^y for Model 1 with Yf 2 = Yf 3 = 0. 
Panel (i) corresponds to benchmark point A, panel (ii) is for B, panel (iii) is for C and 
panel (iv) is for D. The 9 assignments are shown with the separate lines: C 5lS2 (solid), C^ 1 
(dashed), and C^ 1 (dot-dash). The solid curve corresponds to zero D-terms, and the other 
curves correspond to different models for the D-terms. The 2002 experimental limit is 
also given by the horizontal line. 
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Figure 3: Plots showing the Branching Ratio for r — ► ^7 for Model 1 with = Yf 3 = 0. 
Panel (i) corresponds to benchmark point A, panel (ii) is for B, panel (iii) is for C and 
panel (iv) is for D. The assignments are shown with the separate lines: C 5lS2 (solid), Cf 1 
(dashed), and C^ 1 (dot-dash). The solid curve corresponds to zero D-terms, and the other 
curves correspond to different models for the D-terms. The 2005 experimental limit [27] is 
also given by the horizontal line. 
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Figure El shows analogous results for BR(r — > /ry) for Model 1, plotted against 
the gravitino mass 77^3/2- All benchmark points come below the experimental limit 
for a substantial amount of the parameter space. The experimental limit here is more 
recent, and subsequently much more stringent than the previous limit. For these 
models in which there is a large (2,3) element in the neutrino Yukawa matrix the 
branching ratio for r — ► /17 is essentially as constraining as that for /1 — > &y, as 
first pointed out in [§]. The D-term coupling to right-handed scalars has a Yukawa 
mixing angle of order A 3 , compared to A 2 for \i — > &y. A « 0.22 is the Wolfenstein 
parameter, which contributes on an equal footing to e and 5. So the right-handed 
sector is of equal importance to the left-handed sector. We note that the see-saw 
effect enters prominently in the left-handed sector, and by considering Eqs. (|51|) and 
(|52|) for the soft scalar masses in the (2 — 3) sector for r — > fi-f, one can show that 
there is little effect coming from the D-terms coupling to left-handed scalars, since 
we have universal family charges for the left-handed (2 — 3) sector, q^ 2 = qi 3 . For the 
right-handed scalars, however, the D-terms do play an important part and can have 
rather interesting and surprising effects, as we now discuss in some detail. 

The solid line in panel (ii) of Figure H3 for the C 5l5i string assignment of 9 has 
zero contribution from the U(l)p D-terms, and shows just the effect of non-minimal 
SUGRA. This actually suppresses the flavour violation arising from the see-saw effect 
alone, showing an interesting cancellation between the LFV from the see-saw mech- 
anism and the LFV from the non-universal D-terms. On the other hand the dashed 
line in panel (ii) for the C\ 1 case is very similar to the see-saw scenario of benchmark 
point A shown in panel (i). This is due to the D-terms actually conspiring to restore 
universality in the scalar masses, turning non-minimal SUGRA back into the minimal 
form. One can easily see this by applying Eq. (|fi2j) to Eqs. (|53*|) - f)5fi|) for the right- 
handed scalar mass matrices, as the D-terms bring in an equal but opposite effect to 
the non-universal effects from SUGRA, and subsequently force the mass matrices to 
become universal. It is an amazing consequence of this string assignment for 9 that 
in this model the effects of the non-universal U(l)p D-terms can exactly cancel the 
effects of the non- universal SUGRA for the branching ratio of r — > /17, leading to 
universal scalar mass matrices, even with SUGRA turned on. This is a string effect 
that directly affects the amount of flavour violation predicted in this scenario. For 
the Cl 1 case shown by the dot-dash line, applying Eq. ()63|) to the right-handed scalar 
mass matrices Eqs. f|53[) - ()56|) shows that the D-terms in this case actually enhance 
the effect of non-minimal SUGRA, causing the scalar mass matrices to become even 
more non-universal. One can understand this purely right-handed effect by consider- 
ing the different mass insertion diagrams for the left- and right-handed sectors. The 
left-handed sector involves charginos, whereas the right-handed sector only involves 
neutralinos, so the right-handed masses scale differently with the gravitino mass as 
compared to the left-handed masses. Thus we do not have a universal mass scaling 
between the left- and right-handed sectors. This leads to the observed smooth cancel- 
lation of flavour violation between the two competing contributions from the see-saw 
mechanism and from SUGRA with additional U(1)f D-terms. 

the corrected values of S, T3, as discussed. 
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Figure 4: Plots showing the Branching Ratio for fi — > for benchmark point B for Model 
1 only. Panel (i) has Yfe = and Yf 3 = 0. Panel (ii) has Yfc = 1.5 x 10" 3 and Yf 3 = 0. Panel 
(iii) has = and = 1.5 x 10~ 2 . Panel (iv) has Yfa = 1.5 x 1(T 3 and Yf 3 = 1.5 x 10~ 2 . 
The assignments are shown with the separate lines: C 5l5i (solid), Cj* 1 (dashed), and C^ 1 
(dot-dash). The solid curve corresponds to zero D-terms, and the other curves correspond 
to different models for the D-terms. The 2002 experimental limit |26j is also given by the 
horizontal line. 



Figure |3] shows benchmark point B for Model 1. The four panels show the Yf 2 and 
Y± 3 electron Yukawa elements being turned on and off. The results for Figure Efi) 
are for Yf 2 = and Y{ 3 = 0. This is the same as in panel (ii) of Figure El and is the 
base from which we start. Panel (ii) of Figure 0] has Yf 2 = 1.5 x 10 -3 and Y{ 3 = 0, 
so we can clearly see the effect of turning Y± 2 on. It only affects the C 5l5:2 line, as the 
D-terms dominate over this effect for the other two string assignments. Panel (iii) of 
Figure |U uses Y* 2 = and Y± 3 = 1.5 x 10~ 2 , highlighting the effect of just Y± 3 alone. 
Again the zero D-term line of C 5lS2 is the only one that is sizably affected by this 
change in Yukawa texture. Panel (iv) of Figure 0] shows the effect of turning on both 
Yukawa elements: Yf 2 = 1.5 x 10 -3 and Yf 3 = 1.5 x 10~ 2 . We see that the shape of 
the solid line is determined by both Yukawa textures - they seem to have an equal 
impact on it. 

Figure shows benchmark point E, which combines the features of benchmark 
points B and C, thus both U(l) D-term and Froggatt-Nielsen (FN) flavour violation 
appear in the predicted branching ratios shown in this figure. There is some interest- 
ing interplay between the FN F-terms in the benchmark point C region of parameter 
space, and the D-terms in benchmark point B, and benchmark point E is designed 
to show this difference. In Figure the results are shown for Model 1 with both Y± 2 
and Yy 3 Yukawa elements on and off. The shape of the curves are as in benchmark 
point C due to Xg being turned on, and the results are numerically similar to those 
of benchmark point C, showing that the dominant contribution comes from Froggatt- 
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Figure 5: Plots showing the Branching Ratio for \x 



1 only. Panel (i) has Y{ 2 = and Yf 3 = 0. Panel (ii) has = 1.5 x 10~ 3 and Y,% 



ej for benchmark point E for Model 
12 — u aiiu x 13 — u. x auoi ynj naS = 1.5 x 10 -3 and Yf 3 = 0. Panel 
(iii) has Y& = and = 1.5 x 10~ 2 . Panel (iv) has Yf 2 = 1.5 x 1(T 3 and Yj| = 1.5 x 10~ 2 . 
The 6 assignments are shown with the separate lines: C 5l5i (solid), Cj* 1 (dashed), and C^ 1 
(dot-dash). The solid curve corresponds to zero D-terms, and the other curves correspond 
to different models for the D-terms. The 2002 experimental limit |26j is also given by the 
horizontal line. 



Nielsen flavour violation. However, as in benchmark point B, the different D-terms 
corresponding to different 9 assignments leads to noticeable shifts in the results. In 
panel (i) with Y± 2 = and Yf 3 = it is seen that the presence of non-zero D-terms 
actually reduces the LFV rate somewhat compared to the solid curve with zero D- 
terms, corresponding to a region of parameter space where there is some cancellation 
between the flavour violation from the Froggat-Nielsen fields and that caused by the 
U(l)p D-terms. The other panels show variation of the Yukawa elements as in Fig- 
ure EJ for example panel (iv) corresponds to both Y^ 2 and Yf 3 Yukawa elements being 
non-zero. Note that the solid curve in panel (iv) is slightly lower than the solid curve 
in panel (ii), showing that sometimes a non-zero Yukawa coupling can reduce LFV. 

Figure El shows the effects of the Yukawa elements for benchmark point E using 
Model 2, where Model 2 has the same Yukawa structure as Model 1, but has the 
feature that the left-handed family charges are the same for all three families, resulting 
in universal D-terms, at least in the left-handed sector. Comparing Figure El to 
Figure we see that in panel (i) of both figures with Y{ 2 = and Yf 3 = there is no 
obserable difference between the predictions of the two models. However comparing 
panels (iv) of Figure Inland Figure El we see that with non-zero Yf 2 and Yf 3 Model 2 
has the effect of reducing the LFV resulting from the D-terms. 

We could have presented similar plots for benchmark point D, with a new point 
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2 only. Panel (i) has Y{ 2 = and Yf 3 = 0. Panel (ii) has Yf 2 = 1.5 x 10~ 3 and Yf, 



ej for benchmark point E for Model 
12 — u aiiu x 13 — u. x anci V n; naS = 1.5 x 10 _3 and Yf 3 = 0. Panel 
(iii) has Y& = and = 1.5 x 10^ 2 . Panel (iv) has Yf 2 = 1.5 x 1(T 3 and = 1.5 x 1(T 2 . 
The 6 assignments are shown with the separate lines: C 5l5i (solid), Cj* 1 (dashed), and C^ 1 
(dot-dash). The solid curve corresponds to zero D-terms, and the other curves correspond 
to different models for the D-terms. The 2002 experimental limit |26j is also given by the 
horizontal line. 



F, 16 which is an amalgamation of points B and D, but since this produces the same 
kind of results as Figure El we did not present the results. 



7 Conclusions 

We have catalogued and quantitatively studied the importance of all the different 
sources of LFV present in a general non-minimal SUGRA framework, including the 
effects of gauged family symmetry. We have discussed five different sources of LFV in 
such models: see-saw induced LFV arising from the running effects of right-handed 
neutrino fields; supergravity induced LFV due to the non-universal structure of the 
supergravity model; FN (Higgs) flavour violation, due to the F-terms associated with 
FN (Higgs) fields developing vevs, and contributing in a non-universal way to the 
soft trilinear terms; D-term flavour violation, where the D-term mass correction from 
the breaking of the Abelian family symmetry drives the scalar mass matrices to be 
non-universal; and finally the effects of different choices of Yukawa textures on LFV. 

16 The values for the X parameters for benchmark point F would be 

Xs Xt Xh Xjj Xg Xg 

0.290 0.264 0.595 0.595 0.000 0.000 ' 
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In order to quantify the importance of the different effects we investigated these 
disparate souces of LFV numerically, both in isolation and in association with one 
another, within a particular SUGRA model based on a type I string-inspired Pati- 
Salam model with an Abelian family symmetry, which has a sufficiently rich structure 
to enable all of the effects to studied within a single framework. Within this framework 
we derived the soft supersymmetry breaking terms, including the effect of the D- 
terms associated with breaking the family symmetry. For these models the D-terms 
are calculable, but are model dependent, depending on a particular choice of string 
assignment for the FN fields, and in particular the D-terms are only non-zero for 
the non-universal SUGRA models. We have performed a detailed numerical analysis 
of the five sources of LFV using five benchmark points designed to highlight the 
particular effects, and we have explored the effect of the variation of Yukawa texture 
elements on the results. 

The most striking conclusion is how dangerously large the calculable D-term con- 
tribution to flavour violation can be, at least for the class of models studied. However 
it should be emphasised that while the D-terms are calculable in these models they 
are also model dependent, and it is always possible to simply switch off the D-terms 
by selecting the 9 field to have the same string assignment as the 9 intersection state. 
However other choices will lead to non-zero but calculable D-terms, which can be dan- 
gerously large, or can massively suppress flavour violation. For example the curves 
with non-zero D-terms in Figure 0] all exceed the experimental limit for the Branching 
Ratio for fi — > ej, showing that D-term effects have the potential to greatly exceed the 
other contributions from SUGRA, the see-saw mechanism, FN and Higgs, depending 
on the choice of Yukawa textures. However in some cases the D-terms generated 
by breaking the U(1)f family symmetry can also suppress the Branching Ratio for 
fi — > ey, as shown in panels (i) and (iii) of Figures El and H3 

Another notable feature is the effect of Yukawa texture on the results. The Yukawa 
texture has Yf 2 = Y^ 3 = 0, and we have shown that turning on non-zero values of 
these Yukawa couplings can greatly enhance the branching ratio for /i — > almost 
arbitarily. The reason is that the rotations to the SCKM basis are controlled by these 
Yukawa elements and the larger these rotations the larger will be the off-diagonal soft 
masses in the SCKM basis. The non-zero magnitudes of Y{ 2 and Y{ 3 were therefore 
chosen to be large enough to show the variations in the branching ratios of the different 
models, but small enough to keep within the currently experimentally allowed range. 

In this paper we have worked with a particular Yukawa texture in which there is 
a large (2,3) element in the neutrino Yukawa matrix leading to large see-saw induced 
LFV and a branching ratio for r — > fi'-f which is as constraining as that for fi — > ej 
However we have seen that in some cases the D-terms can lead to a large suppression 
of the rate for particular values of m 3 / 2 , as seen in panel (ii) of Figure El For other 
cases the effects of the non-universal U(l)p D-terms can exactly cancel the effects of 
the non-universal SUGRA model leading to universal scalar mass matrices, thereby 
restoring universality even for a non-minimal SUGRA model. Such effects are only 
possible in certain string set-ups and thus LFV provides an observable signal which 
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may discriminate between different underlying string models. 

In conclusion, we have seen that within realistic non-minimal supergravity models 
there can be several important effects leading to much larger LFV than in the case 
usually considered in the literature of minimum flavour violation corresponding to 
just mSUGRA and the see-saw mechanism, and considered here as benchmark point 
A. We find that the D-term contributions are generally dangerously large, but in 
certain cases such contributions can lead to a dramatic suppression of LFV rates, 
for example by cancelling the effect of the see-saw induced LFV in r — > /ry models 
with lop-sided textures. In the class of string models considered here we find the 
surprising result that the D-terms can sometimes serve to restore universality in the 
effective non-minimal supergravity theory. Thus D-terms can give very large and very 
surprising effects in LFV processes. In general there will be a panoply of different 
sources of LFV in realistic non-minimal SUGRA models, and we have explored the 
relative importance of some of them within a particular framework. The results here 
only serve to heighten the expectation that LFV processes such as fi — » and 
t — > ^7 may be observed soon, although it is clear from our results that the precise 
theoretical interpretation of such signals will be more non-trivial than is apparent 
from many previous studies in the literature. 
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A Parameterised trilinears for the 42241 model 

We here write the general form of the trilinear parameters assuming nothing 
about the form of the Yukawa matricies. 

\/3m 3/2 {X s [l + {S + S) d s In Y abc ] 

+ X Tl [-l + fr + Tj^hiYate] 
+ X T2 [-1 + (T 2 + T 2 ) d T2 \nY a J] 
+ X T3 (T 3 + T 3 )d T:i \nY abc 

+ X H (S + S)^Hd H lnY abc 

+ X 7T (T 3 + T 3 yHd w \nY abc 

+ X e (S + Sy (T 3 + T 3 ) 1 6d e In F a6c } (A.l) 



^c 5l c 5l c 5 i 5 2 = V / 3m 3 / 2 s X, 

1 3 



-+ (S + S)d s In Y abc 

X Tl [-1 + (Tx + Tx)^ In y afec ] 
X T2 (T 2 + T 2 )d T2 \nY abc 



+ X_ 



T 3 



~(T 3 + T 3 )d T3 lnY abc 



+ X H (S+S)*Hd H \nY abc 

+ X 7T (T 3 + T 3 yHd w \nY abc 

+ X e (S + S)^ (T 3 + T 3 ) * 9d e In F abc } 



(A.2) 



~ + (S + S)d s lnY abc 

+ X Tl [-1 + (T 1 + T 1 )a ri In Y abc ] 
+ X T2 [-1 + (T 2 + T 2 ) 9 T2 In F afec ] 



X. 



T 3 



-(T 3 + T 3 )d T:j lnY abc 



+ X H (S + Hd H \nY abc 

+ X 7T {T 3 + T 3 YHdj I \nY abc 

+ X e (S + Sy (T 3 + T 3 ) 1 9d e In F abc } 



(A.3) 
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A c^c^c Sl = ^ m 3/2 {X s (S + S) d s In Y abc 

+ X Tl [-l+^+f^lnY^J 
+ X T2 (T 2 + T 2 )d T2 \nY abc 
+ X n (T 3 + T 3 )d T3 \nY abc 

+ X H (S + S) h Hd H \nY abc 

+ X 7T (T 3 + T 3 yHd w lnY abc 

+ X e {S + Sy (T 3 + T 3 ) 1 6d e In F abc } (A.4) 

B Derivation of the D-terms 

What follows is a full derivation of D H and D e from the superpotential of the U{l)p 
extended supersymmetric Pati-Salam 42241 model, as detailed in Section 0] 

The relevent parts of the superpotential for the 42241 model are those concerning 
the Higgs and Froggatt-Nielsen fields which have different gauge singlets 17 . 

W = SX S (HH - M 2 H ) + S'\ S '(60 - M 2 ), (B.l) 

where the Mjj and Mg are GUT scale masses associated with the Higgs and Froggatt- 
Nielsen vevs respectively. We have assumed that the heavy Higgs develop vevs along 
the neutrino directions only, such that 

(H) = H U = H ; (H) = H U = H. (B.2) 

This is because charged objects gaining vevs would break their charge group at 
the GUT scale, causing problems 18 which the neutral components avoid. Similarly, 
the Froggatt-Nielsen vevs are concisely written as 

(6) = 6 ; (6) = 6. (B.3) 

The F-terms associated with the singlet fields are 

\F S \ 2 = \X S (HH - M 2 H )\ 2 (B.4) 

and 

\Fsj 2 = \M99 - M 2 )\ 2 . (B.5) 

17 Notc that these must still have the same quantum numbers as they are both singlets, and 
therefore in the same representation of the gauge group. 

18 A colour charged object, for example, would imply that QCD is broken at the GUT scale, which 
would lead to a very massive gluon. 
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We use these to form our Higgs potential, 

V H = V D + V F + V soft . (B.6) 



The F-term potential is trivially obtained from Eqs. (jB.4|) and (jB.5|) . and the soft 
terms are simply written down with mass-squared terms for each of the soft SUSY- 
breaking scalar masses associated with the Higgs and FN vevs, as we will later see. 
The D-term potential takes a little more work, so we shall cover that here. The 
general form for the D-term potential is 

1 1 3 1 15 

Vd = z&DIDf + -/2R E WrWr + 2 91 £ (B.7) 

a=l m=l 

where gp is the gauge coupling for U{1) F , 9ir is for SU(2)r, and is for 577 (4) c . 
We focus on the a = 3 and m — 15 contributions to Eq. (jB.7j) which involve the 



and 



■it = dia S I 5. -5 I (B.8) 



the diagonal generators of the SU(2)r and S77(4) groups. There is no sum in the 
£7(1) part of Eq. (|B.7j) because this group is the unit matrix 1, so the only generator 
is 

T F = lq F . (B.10) 

where q F is the charge of the U(1)f group for each field. When applying Eqs. (JB.8j) 
and ()B.9j) to conjugate fields we have to complex conjugate the generator and multiply 
it by —1, but for the U(1)f group we just use Eq. (jB.lOJ) where it is known that the 
qF charges are different for right-handed fields. 

Using Table H] we find that Dp, D^, and D\ 5 are given by 

Dp = H\q F )HH + H\q F ) H H +tf(q F )rf + 6\q F ) e 

+ F\{q F )pF 3 51 + F}{qp) Fi F^l + h\q F ) h h (B.ll) 

Dl R = H\-r^)H + H\rl)H + F\-T^)F + h\Tl)h (B.12) 

D\ 5 = H\-Tl 5 *)H + H\Tl b )H + F j {-Tl 5 *)F + F f (T 4 15 )F. (B.13) 

In Eq. (jB.llJ) . i,j G {1, 2, 3} are family indices. We used Si to pick out the trace 
of the outer products F-F^ and f\f\ thereby giving us the dot product. 

The scalar components of the left-handed matter superfield F are q and /. The 
scalar components of the right-handed matter superfield F are u c , d c , u c , and e c . 
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The tensorial conventions are shown in [22] for D\ R and D\ 5 . For Dp, all diagonal 
elements are equal to unity for the generator, so the tensor notation is trivial. 

Now we square Dp, -D| R , and D\ 5 (as they are squared in the D-term potential), 
then consider the cross terms, as we are only interested in those terms that look 
like (mass) 2 ■ (field) 2 . The (mass) 2 terms come from the vevs of the heavy Higgs 

2 2 

and Froggatt-Nielsen fields H ,H 2 ,6 ,6 2 . The (field) 2 terms come from \Fi\ 2 , \Fi\ 2 , 
where the % is a family index running from 1 to 3. 

In Eqs. (JEU) to (jHHjl we have defined D 2 H and Dj to be 

D\ = \(H 2 -H 2 ) (B.14) 
D 2 e = -q e (f-9 2 )-q H (H 2 -H 2 ), (B.15) 



where qg is defined to be —1 and qu is defined to be —qm, thus qu = \ for Model 
1 and qu = 1 for Model 2. 

By considering the charge structure in Table El we can clearly see how the charge 
structure effects the additional U(1)f D-term contributions to the sparticle mass 
matrices in Eqs. (J5T|) - (J58I) . 

We have followed j22] in choosing our designation of D 2 H , but D 2 e is new. Working 
explicitly within Model 1, where qu = |, we shall proceed to rewrite these D-terms 
as functions of the soft SUSY breaking masses and gauge couplings. So we arrive at 
our D-term potential, which, when put together with Vp and V^ Q f t , forms our Higgs 
potential: 



T/ 1 2 

^Higgs — 2 9f 



2 1(H 2 _ E 2 ) 2 + (f - 9 2 ) 2 - h -(H 2 - H 2 )(f - e 2 



+l(9ln + \gl)(H 2 - H 2 ) 2 + \ 2 S (HH - M 2 H ) 2 + \ 2 S ,(66 - M 2 ) 2 



+m^H z + m 2 H H 2 + m 2 ^ + m 2 6 2 , (B.16) 

where Vd are the terms multiplied by gauge coupling factors, Vp are the terms mul- 
tiplied by the dilaton lambdas, and V^ ft are the last four terms multiplying the TeV 
scale soft SUSY breaking scalar masses. Mu and Mg are GUT scale masses. 

To find the form of the D-terms, we must minimise this potential with respect to 
the fields H, H, 9, 9, and then set these minimisation relations equal to zero. 

As these are set to zero, any linear combination of them is also zero, so taking the 
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following combinations and rearranging them, we have two minimisation conditions 



OV dV 
dV dV 



25 2 i 2 ,^2 

"g- 9f + 92R + 



1 
1 

\g 2 F {e 



(H + Ey - y s {HH - M Z H ) \{H-H) 



6){H + H) = -mLH + m 2 H H (B.17) 



- 6 9l 

-A 5 , 



6 



- 6 2 ) - 5{H - H 2 ) 
Mjj)(6 -6) = -mW + m%6. 



:B.18] 



2 2 

For D-flatness, it is necessary to set H = H 2 and 8 = 9 2 , which results in 
Vd = 0. For i^-flatness and i^z-flatness, it is necessary to set HH = Mjj and 
69 = Mg, yielding zero valued F-terms. We wish to perturb away from these flatness 
conditions, so we impose 



H = M H -m ; H = M H - m 
6 = M e -m ; 6 = M e -m', 

where m, m, m', and m! are all TeV scale masses. 

Thus the two minimisation conditions, Eqs. (jB.17|) and (jB.18|) . become 19 

10 



(B.19) 
(B.20) 



g 2 H (m -m)M H - -jg 2 F {rn' - m')M e 
2 

-g 2 F [6(m' - m')M e - 5(m - m)M H ] 
3 



m 



H 



m — rn-pr 



;b.2i) 
;b.22) 



Now putting the small perturbations, Eqs. (|B.19|) and ()B.20|) . into Eqs. (jB.14|) 
and (|B.15|) for the D-terms, we have, to leading order, 



n 2 



Da 



-(m — ra)Mu 

- [6(m' - m!)M e - 5(m - m)M H } . 
3 



(B.23) 
(B.24) 



So, using; Eqs. l|BT2T) and (jB~22|l in the above Eqs. (jR23jl and (|B~24|l . we have 
the following expressions for our D-terms for Model 1, as functions of the soft SUSY 
breaking masses, GUT scale masses and gauge couplings 



n 2 



1 



±92 R + m 

2 2 

mi — m4 



2 2 , u / 2 2n 



2^ 



(B.25) 
(B.26) 



Note that Eq. (|B.26|) was used in obtaining Eq. (|B.25|) . This is the form of the D- 
terms as used in the updated version of SOFTSUSY j2H] to compute the slepton mass 



After rearranging and taking the leading order in the GUT scale masses Mh and Mg 



35 



data for our lepton flavour violating branching ratios in Model 1, using Eqs. (|51|) to 

dHHD- 

For Model 2, the derivation is very similar, with the factor of qu being the only 
difference, and so we obtain the form of D 2 H below, with D 2 , being the same in both 
models. 



(m 2 M — mZj) + (mi — m|) 

D " = H^A • (&2T) 

In both cases, we can see that the Pati-Salam limit is obtained when the Froggatt- 
Nielsen scalar masses are degenerate, m 2 g = m|. This result differs from [23] due to a 
different derivation procedure. 



C Operators for the two models 

In this appendix we give the operators which are responsible for generating the 
Yukawa matrices of Models 1 and 2 analysed in this paper. 



C.l Model 1 

This model is almost the model studied in |2*Ul Hlj . but with an extra operator in 
the (1,2) and (1,3) Yukawa matrix elements to allow a non-zero Yf 2 and Yf 3 . The 
operator texture is 



O = 



( ail O Fc + a l { 1 O llAe )e> (a 12 Ee + a' 12 0' cb + a'( 2 0" Ec )e 3 (a u O Ec + a' 13 0' c f + a'{ 3 0" Ee )e 
(a 21 Dc )t A (a 22 Bc + a' 22 0' F f)e 2 (a 23 Ee + a' 23 0' Bc ) 

(a 3 i0 Fc )e 4 (a 32 Ac + a' 23 O lFe )e 2 a 33 

(C.l) 
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Model 1 


Model 2 


X 



u.zz 


X 




n 99 
v.aa 


6 


n 99 
u.zz 


6 


n 99 
v.aa 


a 33 




a 33 


u.00 


an 


n n9 


„/ 

a n 


n no 
-U.yz 




U.00 


„/ 

a 12 


U.00 


Ol3 


n nn 
u.uu 


^13 


n nn 
u.uu 


^21 


1 fi7 
1 .0 ( 


«21 


1 fi7 
1 .0 ( 


«22 


119 
1 . LA 


«22 


119 
1 . LA 


u 23 


8Q 


u 23 


89 


«31 


-0.21 


«31 


-0.21 


a 3 2 


2.08 


«32 


2.08 


<2 


0.77 


<2 


0.77 


«'l3 


0.53 


<3 


0.53 


a 22 


0.66 


a 22 


0.66 


fl 23 


0.40 


a 23 


0.40 


a 32 


1.80 


a 32 


1.80 


<1 


0.278 


a'" 
a ll 


0.278 


<2 


0.000 


"12 


0.000 


<3 


0.000 


/// 
a 13 


0.000 


An 


0.94 


An 


0.94 


Al2 


0.48 


An 


0.48 


Al3 


2.10 


A 13 


2.10 


A 22 


0.52 


A22 


0.52 


A 23 


1.29 


A 23 


1.29 


A 3 3 


1.88 


A33 


1.88 



Table 4: The 0(1) coefficients in Model 1 and Model 2. The values of 013, a'{ 2 (in 
Model 1) and a" 2 ( in Model 2) will be varied, as discussed in the text. 
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where the operator nomenclature is defined in Appendix [D] and Appendix |Ej This 
leads to the following Yukawa textures: 



Y"(M 



x, 



Y d (M_ 



x> 



x, 






. a 3 i|5e 4 
-a 31 Ue 4 



'12 



U 22 r ^° fc 



'■5VE 
a' ^S 2 e 2 

a 32 5° fc 



«33 



-a 



12 



4?5 2 e 



l«22 



"13 V5 



2 5 + a' ^ 



a 32 \/25e 



22 5^' 
2 



'23 



|5 2 

5 



a 3 3 



-Q223 



a'U25 2 e 3 

+ outfit)** 



12 

2 /I 
5 V 5 

a 3 2V / 25e i 



a 13 25 2 e 



'23 

a 33 



|5 2 

5 



Y V {M X ) 



nV25 





3^5 



a 12 25e 3 



2^2 



22 5 v^° 6 
/ 6*2 2 

32 5 ° 6 



a 



For both models we define e and 5 as: 

(0) 



M 



x 



a'( 3 26 3 e 
a 23 2S 

a 3 3 



( mm 

V Ml 



(C2) 



(C.3) 



(C.4) 



(C.5) 



• 



(C.6) 



We take 5 = e = 0.22. 



The values of the arbritary couplin 
values for the Yukawa elements which 
values of Y* 2 and Y* 3 inserted instead 



gs are laid out in Table EJ This gives numerical 
can be used in either model, with the relevant 
of the texture zeros: 



(C.7) 
(C.8) 
(C.9) 
(CIO) 





' 2.159 x 10" 


-06 


5.606 x 10- 


-04 


5.090 x 10- 


03 " 


Y U (M X ) = 


0.000 




1.105 x 10- 


-03 


0.000 






_ 0.000 




6.733 x 10- 


-03 


5.841 x 10- 


01 




-1.661 x 10' 


-04 


-5.606 x 10 


-04 


1.018 x 10 


-02 _ 


Y d (M x ) = 


7.683 x 10- 


-04 


-5.343 x 10 


-03 


1.216 x 10 


-02 




-1.769 x 10' 


-04 


3.133 x 10" 


-02 


3.933 x 10 


-01 




-1.246 x 10' 


-04 


0.000 




0.000 




Y e (M x ) = 


1.537 x 10" 


-03 


2.432 x 10' 


-02 


-3.649 x 10 


-02 




-1.327 x 10' 


-04 


3.133 x 10' 


-02 


5.469 x 10' 


-01 




" 2.159 x 10- 


-06 


1.525 x 10- 


-03 


0.000 




Y V {M X ) = 


0.000 




8.290 x 10- 


-04 


3.923 x 10- 


01 




0.000 




5.050 x 10- 


-03 


5.469 x 10- 


01 



The RH Majorana neutrino mass matrix is: 
M RR {M X ) 



Mi 



33 



A u 5e 8 A 12 5e 6 A 13 Se 4 
A 12 6e 6 A 22 5e* A 23 5e 2 
A 13 5e* A 23 6e 2 A 33 



(en; 
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The numerical values for the Majorana mass matrix are 
M RR (M X ) 



M: 



33 



3.508 x 10 s 
3.686 x 10 9 
3.345 x 10 11 



3.686 x 10 9 
8.313 x 10 10 
5.886 x 10 12 



3.345 x 10 11 
5.886 x 10 12 
5.795 x 10 14 



(C.12) 



C.2 Model 2 

The operator texture for Model 2 is: 



O 



[a' lx O lFc 

(a 2 iO Dc y 
(a 31 Fc )e 4 



a 'l> i0 "'Aey ( a > 12 o>Ee + a '{ 2 0" cb + a'{ 2 0"' Ec )t 2 (a 13 Ec + o![ z O" c ^ + a'{ 3 0"' Ee ) 



(a 22 Bc + a' 22 0' F f)e 2 



(a 32 Ac + a> 23 0' Fe )e 2 



(a 23 Ee + a' 23 0' Bc ) 



«33 



(C.13) 



The operator nomenclature is defined in Appendix [D] and Appendix |E| The new 
operator setup leads to the following Yukawa textures: 



Y U {M_ 



x, 



Y d (M 



x, 



Y e (M_ 



x, 



Y V {M^ 



x , 






~5 2 e A 

2 xA 
i 

6 ^2 A 

4 5e 4 



H5 l 



a21 7E Se 



a n f5 2 e 



021 



4 
4^4 



-a 3 i|5e 






3^2 



12 v ^ c a'{ 3 ^rd 

/ 8 5 2 e 2 



a 



32 5 l 



V 3 

o' 

a 3 3 



-a / ' v / 25 3 e 2 



«22 



'12 



a 32 \/25e 



22 5^5 
2 



'2:; 



U 2 

5 



^33 



'12 



2<5 3 e 



a 13 2<5 3 



'— a 2 23 



'l + a' 22 5^)6e 2 -a> 2? JU 2 



a 32 V25e 2 



'23 
^33 



a' l2 2S 2 e 2 



2 r 2 



a'JPe 2 



a 23 2S 
a 3 3 



(C.14) 



(C.15) 



(C.16) 



(C.17) 



The RH Majorana neutrino mass matrix is the same as in Model 1: 

M RR (M X ) 



33 



A u 5e 8 A 12 5e 6 A 13 Se 4 
A 12 5e 6 A 22 5e* A 23 5e 2 
A 13 5e* A 23 5e 2 A 33 



(C.18) 



D n = 1 operators 



The n = 1 Dirac operators are the complete set of all operators that can be con- 
structed from the quintilinear FFhHH by all possible group theoretical contractions 
of the indicies in 

= F aa F Px hlH lz H? w (D. 1) 
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Operator Name 



-qZZ- 

Ab 

qAc 

QAd 

Q Ae 

QBa 
QBb 

QBc 

QBd 

QBe 

QBf 
QCa 
QCb 
QCc 

QCd 

QCe 
QCf 

QDa 
QDb 

QDc 

QDd 

QDe 

QDf 
QEa 
QEb 

qEc 

QEd 
QEe 
Ef 
QFa 
QEb 

QEc 

QEd 
QEe 
QFf 



Operator Name in j2I] 



B 

QM 

T 

o v 
o u 
o c 

D 

QW 

o x 
o z 

Y 

Qa 

F 

QE 

o b 

QN 

o c 

Qd 

o e 

QG 

H 
0° 

of 

09 
Qh 

l 

oj 
o 1 

J 

o p 

Q 

QR 

L 
K 
s 



QUh 2 



QDh^ 



LEh t 





2V2 
5 

yl 
5 
1 

VE 
1 

V5 



2\/2 

5 

r 2 



5 

>/2 
V2 



4 

2 

2 

VL 
2 
5 

2 
5 



j_ 

5 
2 

V5 
2 
5 








Ay/2 

5 

Ay/2 
5 

16 



-1 
V2 

5 



2V2 

5 
1 

V5 
-1 

V5 

t~2 

5 

5 



2V2 
5 

-V2 

2 

2 



4 



2 

2 
5 



5 








Ay/2 

5 

~AV2 



5 

8 



16 
5^ 



1 

-1 

V2 

5 



2V2 

5 

-3 

-3 
V5 



-3V2 
5 



-6^ 
5 










V 5 



4 

J_ 

5 



8 
5 

V2 

2 

2 



4 

3v^2 
5 

-3>/2 



5^5 



12 

5VE 



Table 5: Operator names, CGCs and names in [21] 
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We define some SU(4) invariant tensors C and some SU(2) invariant tensors R as 
follows 20 : 



(Ci)? 


- °/3 




_ ,P7^X 


l°ioJ P7 


= s a S p + S a S p 

P 7 7 P 


(^15) ctp 






= ^ x 


(As)™ 





Then the six independent SU(4) structures are: 
A- (COJCGOJ = ^ 



B- (C 15 )S(^i 5 )^ = 

C (C 6 )2(C 6 )S = 8(^-5^) 

D. (C 10 )^ (Ci„)£ = 2(^ + 5^) 

e. (co^coi = ^ 

F- (Ci 5 )^(C 15 )^ = (D.3) 



And the six SU(2) structures are: 

b- (Rs)Z(Rs)Z = «-k 



3/9 w y 2 V w 

= e xz e 



C. 6 6 6 ? 

d. e ws e°* (R 3 ) s y l (R 3 y/ q = 6*5* - \e wy ^ 
(R 1 ) z y (R 1 )l = fi*% 



All possible n = 1 operators were then named A ...O z O a ...O j in We rename 
them here in a manner consistent with the n > 1 operators O*™'- 1 , so that the names 
are O nn where II is the SU(4) structure and n is the SU2 structure. See Table El for 
the translation into the names of Ref.j2I] and the CGCs. 

20 Thc subscript denotes the dimension of the representation they can create from multiplying 4 
or 4 with 4 or 4. For example (Cis)^4 7 4 p = 15^ . 
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All of these operators are operators for the case without a U(l) family symmetry. 
In the case when there is, we follow the prescription 

/ e Y ij 

■ (D5) 

Where pu = \Xo ZJ \ is the modulus of the charge of the operator. If the charge 
of the operator is negative, then the field 9 should be replaced by the field 9. The 
prescription makes the operator chargeless under the U(1)f while simultaneously not 
changing the dimension. 



E n > 1 operators 

In the case that n > 1, there will be more indicies to contract, which allows more 
representations, and hence more Clebsch coefficients. To generalise the notation, 
it is necessary only to construct the new tensors which create the new structures. 
However, it will always be possible to contract the new indicies between the H and H 
fields to create a singlet HH which has a Clebsch of 1 in each sector u, d, e, v. In this 
case, the first structures are the same as the old structures, but with extra 5 symbols 
which construct the HH singlet. 

Thus taking an n = 2 operator, say 0' Fb , which forms a representation that could 
have been attained by a n = 1 operator, the Clebsch coefficients are the same. This 
is what we mean by £) n/ri7r , as we have only used n > 1 coefficients which are in the 
subset that have n — 1 analogues. 
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